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Abstract 

The report of a significant deviation of the CMB temperature anisotropics distribution from Gaussianity (soon 
after the public release of the WMAP data in 2003) has become one of the most solid WMAP anomalies. This 
detection grounds on an excess of the kurtosis of the Spherical Mexican Hat Wavelet coefficients at scales of around 
10 degrees. At these scales, a prominent feature — located in the southern Galactic hemisphere — was highlighted 
from the rest of the SMHW coefficients: the Cold Spot. This article presents a comprehensive overview related 
to the study of the Cold Spot, paying attention to the non-Gaussianity detection methods, the morphological 
characteristics of the Cold Spot, and the possible sources studied in the literature to explain its nature. Special 
emphasis is made on the Cold Spot compatibility with a cosmic texture, commenting on future tests that would 
help to give support or discard this hypothesis. 



1 Introduction 

Besides the great success of the NASA WMAP satel- 
lite on providing a detailed knowledge of the cosmo- 
logical parameters that define the physical properties 
of the Universe (e.g., [69, 41 ), some unexpected re- 
sults have attracted the attention of the cosmological 
community soon after the first release of the WMAP 
data: the so-called WMAP anomalies. Some of these 
anomalies are related to hemispherical asymmetries 
(e.g., [la ISg Ha [591 [Ml ESI Ell [ESI [89]), an anoma- 
lous alignment of the quadrupole and octopole compo- 
nents (e.g., [ini [Ml [13 [33 [ni [Til [21 [11 [S^ sig- 
nificantly low variance of the CMB temperature fluc- 
tuations (e.g., [821 HI im [33]) 5 OT anomalous align- 
ment of the CMB features toward the Ecliptic poles 
(e.g., \127\ 1120] ). Some of these aspects are addressed 
in this special issue. In addition to the previous find- 
ings, the prominent cold spot (hereinafter, the Cold 
Spot) detected in the southern hemisphere by [118 be- 
came one of the most studied anomalies of the WMAP 
data. The Cold Spot was identified after testing that 
the Spherical Mexican Hat Wavelet (SMHW) coeffi- 
cients of the WMAP data presented an excess of kur- 
tosis (at scales of around 10° in the sky), as compared 
to the distribution derived from isotropic and Gaussian 
CMB simulations. This paper presents a complete re- 
view on the detection and characterization of this non- 
standard signature, and a description of the different 
attempts made so far in understanding what could be 
the cause behind such departure from the standard in- 



fiationary paradigm. It is organized as follows: in Sec- 
tion |2j I justify the use of wavelets as a natural tool 
for probing the Gaussianity of the CMB temperature 
fiuct nations. I also present the different statistics ap- 
plied in the wavelet space that had led to point out the 
WMAP data incompatibility with the standard model. 
In section [3] I briefiy describe the morphological char- 
acteristics of the Cold Spot. The important question 
of the actual significance of the detection of the Cold 
Spot, and the aspects associated with a posteriori in- 
terpretations are addressed in Section |4| Some of the 
different sources that have been considered in the liter- 
ature to explain the Cold Spot feature are discussed in 
Section |5j In Section [6) I explain in detail a plausible 
hypothesis to accommodate the existence of the Cold 
Spot together with the standard cosmological model: 
a cosmic texture. In addition, I also describe possible 
follow-up tests that could help to confirm or discard 
such hypothesis. Finally, my conclusions are given in 
Section [3 

2 The non-Gaussianity detection 

The Cold Spot was firstly identified through a blind 
Gaussianity test of the WMAP first year data [118 . 
This test was designed to probe the isotropic and Gaus- 
sian nature of the CMB, as predicted by the standard 
infiationary model (see, for instance [72 ), and it was 
based in a multiresolution analysis performed with the 
Spherical Mexican Hat Wavelet (SMHW). In this Sec- 
tion, I summarize the non-Gaussianity detection that 
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led to the identification of the Cold Spot. I start by jus- 
tifying why an analysis based on wavelets was proposed, 
and presenting the main characteristics of the wavelet 
used in the analysis: the SMHW. Afterwards, I explain 
which statistics (all of them based on the SMHW co- 
efficients) reported the original deviations from Gaus- 
sianity of the WMAP data. 

2.1 Why a wavelet? 

Nowadays, the CMB scientific community is already 
very familiar with the application of wavelets (and other 
members of the -lets zoo, like curvelets, ridglets or 
needlets) to data analysis. However, it is worth re- 
calling that this is a relatively new custom. Although 
wavelet applications in cosmology shyly started already 
in mid 80s, it was not until 1997 that the first appli- 
cation to CMB was presented [46], precisely in a ex- 
ercise devoted to probe the Gaussianity of the CMB; 
and it was just a year afterwards that the first CMB 
data analyses with wavelets came to light [971 EO] ^ in 
particular, with two applications to COBE [8 data. 
During the last decade, the application of wavelets to 
extract cosmological information from CMB data has 
growth considerably, in many different branches: com- 
pact source detection (e.g. [ml [HI [ml [5T | l74| l96]\ 
Gaussianity (e.g. P [84l ^ [T9l [78l ITTSl [28l [7^ [Ml IT^ 
[34t llQlj ). cross-correlation with large scale structure 
(e.g. |119[ l93l [73l [8T]), decomposition of the coupled 
E/B signals [16], probing isotropy (e.g. (127111201 155]). 
cosmic string detection [57 , microwave sky recovery 
(e.g. [75l[83l|60l|36 ), CMB denoising (e.g. [Mllm]), 
CMB power spectrum determination [44], and primor- 
dial power spectrum recovery (e.g. [85]) are some of the 
application fields. I refer the reader to [80 11121] for some 
reviews on the wavelets applications to CMB data, with 
a particular emphasis on data analysis on the sphere. 

The wavelet transform (e.g. [88^ ^) has become 
very popular for a major reason: they offer a unique 
opportunity to probe scale-dependent phenomena, but 
keeping, at the same time, information about spatial 
localization. This is a clear advantage — for many 
purposes — over classical Fourier or harmonic trans- 
forms: physical processes typically exhibit a clear scale- 
dependent behaviour, and, often, such behaviour differs 
enough from one phenomenon to another (e.g., in a mi- 
crowave image, the localized emission due to cluster of 
galaxies has very different properties as compared to 
the large-scale signal produced by the Galactic compo- 
nents) . 

The capability of emphasising or amplifying some 
features (at a particular scale) makes wavelets unique 
to probe the Gaussianity of the CMB: there are many 
different physical processes that might introduce non- 
Gaussian signatures into a CMB signal, at a very partic- 
ular scale range (e.g., primordial non- Gaussianity due 



to non-standard infiationary scenarios, cosmic defects 
like strings or textures, secondary anisotropics, fore- 
ground emissions, . . . ). 




Figure 1: Figure shows the profile on the SMHW as a 
function of the angular distance for different wavelet 
scales ranging from 1° to 45°. Notice that the am- 
plitude of the SMHW at l9 = has been fixed to 1, to 
allow for an easier comparison. 

In this sense, when a Gaussianity analysis of the 
CMB is performed in wavelet space, we are in a very 
adequate framework to probe, almost separately, any 
potential non-Gaussian signatures present in the data. 
The ability of wavelets to amplify a given signature is 
explained because they can be seen as compensated fil- 
ters (i.e., J^^ ^ = 0, where is the space where the 
wavelet ^ is defined). This property is not satisfied 
by other standard filters, as the top-hat or the Gaus- 
sian functions. As mentioned above, wavelets do not 
only provide us with the capability of selecting a given 
scale range, but they also allow one to keep the in- 
formation about spatial localization: we do not just 
study, for instance, the compatibility with Gaussianity, 
but we are also able to identify where in the data a 
deviation might be spatially localized. This intrinsic 
property of the wavelet transform (spatial localization) 
is also unique to explore not only Gaussianity, but also 
isotropy, since the statistical properties can be studied 
(almost independently) from one region in the data to 
another, in a self-consistent way. Hereinafter I will fo- 
cus in the use of the two-dimensional (2D) continuous 
wavelet transform (CWT, see, for instance. Chapter 2 
in [4 ), and, in particular, on the 2D CWT defined on 
the sphere. There are different ways to define wavelets 
on the sphere (e.g. [3l ll26[fTQ5lll28| ): most of the works 
related to the Cold Spot have been performed using 
the definition of the isotropic Spherical Mexican Hat 
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Wavelet (SMHW) proposed by [78| which is a stereo- 
graphic projection of the Mexican Hat Wavelet, as pro- 
posed by |3]: 



2ttRNr 



'-(1)1 [-( 



(1) 

where y = 2tan| is the stereographic projection vari- 
able, G [0, tt) is the co-latitude, and the constant 

Nr : 



^ 1 + ^ + ^ is chosen such as the square of 
the wavelet function {0; R) is normalized to unity. In 
figure [l] the radial profile of the SMHW, for different 
wavelet scales i?, is shown. For a signal T (^, <p) defined 
on the sphere, its spherical harmonic coefficients t£rn 
are defined as: 
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where dfl = d^sin^d^, the spherical coordinates are, 
as mentioned above, the co-latitude (related to the 
latitude 6as6 = 7r/2 — ^), and (j) G [0, 27r) is the longi- 
tude. The function Y^rn is the spherical harmonic 
of order £ and m, and * denotes complex conjugation. 
For an isotropic 2D CWT on the sphere the wavelet 
coefficients w (^, R) are obtained as: 



w (l9, (j); R) 
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where {R) is the window function associated with 
the wavelet function (e.g., ^ {0]R) in equation [T]), and 
^max represents the maximum multipole associated with 
a given resolution of the signal T(^,0), typically lim- 
ited by the size of the pixel adopted to represent such 
signal on the sphere. 



2.2 The statistics 

The Cold Spot was firstly detected via a positive devi- 
ation of the kurtosis of the SMHW coefficients at scales 
of around R ~ 300°. The inspection of the map of the 
SMHW coefficients at these scales revealed the presence 
of a very large and cold spot in the southern hemi- 
sphere. The comparison of the amplitude and the area 
of this cold spot as compared with Gaussian simulations 
showed that it was particularly anomalous. Finally, a 
higher criticism test of the SMHW coefficients also in- 
dicated a deviation, at around the same wavelet scales, 
also showing that the major source for such deviation 
was located in the position of the Cold Spot. In the 
following subsections I describe briefly these statistics: 
the kurtosis, the amplitude, the area and the higher 
criticism. 
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Figure 2: Figure showing the positive deviation of the 
kurtosis of the wavelet coefficients at scales R of 250 and 
300 arcmin found by |118j . Red, green and magenta re- 
gions represent the acceptance intervals at 32%, 5% and 
1%, respectively. 

2.2.1 The kurtosis 

The two most obvious indicators for a possible de- 
viation from Gaussianity of a given data set z = 
{zi, Z2, . . . , ^^at} are the third (/is) and fourth cen- 
tral moments. The n-central moment of a random set 
of numbers is defined as: 



1 ^ 

n = ^ (^i - ^) 7 
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where z = (l/N) YliLi mean value of the data 

set. For a Gaussian random data set of zero mean (i.e., 
/ii = z = 0) and dispersion a (i.e., /i2 = cr^), it is trivial 
to prove that higher order central moments are either 
zero (in the case of /i2n+i, for n > 1) or a given func- 
tion of the dispersion (in the case of /i2n7 foi" ^ ^ !)• 
Usually, it is much more convenient to work with nor- 
malized central moments Un, such as: 



(5) 



The normalized central moments are more convenient 
than central moments, since they are referred to the in- 
trinsic fluctuations of the random data set (represented 
by the dispersion). This helps to absorb into Vn possi- 
ble uncertainties on the knowledge on the amplitude of 
the fluctuations of the random sample. The normalized 
central moment is normally referred to as the skew- 
ness (5), whereas = — 3 is called kurtosis. The 
reason for the subtraction of the number 3 in the previ- 
ous expression comes from the fact that, for a Gaussian 
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Figure 3: Left panel: CMB cleaned map derived from the 4tli release of the WMAP data, obtained via the template 
fitting technique described in [50 . Right panel: wavelet coefficients of the previous map, obtained after the SMHW 
convolution at a scale of = 250 arcmin. The location of the Cold Spot is indicated in both panels by the circle. 
The centre of the Cold Spot is = 147° and ^ = 210°. 



random variable, /i4 = 3cr^ and, therefore, = 3. The 
previous definition assures that, as it happens for the 
skewness (5), the kurtosis {K) of a Gaussian field is 
zero. Applying these concepts to the SMHW coeffi- 
cients w (^, (/); i?), the skewness Sr and the kurtosis Kr 
of the wavelet coefficients, as a function of the wavelet 
scale R, can be defined as follows: 
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where, at each scale, it is assumed that the coefficients 
w {6^(j); R) have zero mean. aR is the dispersion of the 
wavelet coefficients at the scale R: 
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In the previous expressions, Npix(i^) represents the 
number of wavelet coefficients at a given scale R. Notice 
that, very often, CMB data is not available in the full 
celestial sphere (for instance because strong contamina- 
tion from astrophysical foregrounds has to be masked). 
In most of the works in the CMB field, the HEALPix 
tessellation [52 is adopted. In this scheme, the resolu- 
tion of a given image represented on the sphere is given 
by the Nside parameter, which indicates how many di- 
visions of the 12 basic pixels are required to achieve 
such resolution. The Nside parameter is related with the 
number of the pixels (Npix) required to fill the sphere 
at that resolution as: Npix = 12Nside ^. In figure [2]the 
major result of the seminal work [118 on the Cold Spot 
is shown. It represents the kurtosis of the wavelet co- 
efficients, as a function of the scale, for the first release 
of the WMAP data (blue stars). The sohd yellow line 



represents the mean value obtained from 10,000 CMB 
Gaussian random simulations, taking into account the 
instrumental properties of the analyzed data, and gen- 
erated from an angular power spectrum derived form 
the best-fit cosmological model. The coloured regions 
(red, green and magenta) represent the acceptance in- 
tervals at 32%, 5% and 1%, respectively. Notice that 
at R scales of 250 and 300 arcmin, the kurtosis of 
the WMAP data was above the 1% acceptance inter- 
val. In detail, the excess of kurtosis is given by a p- 
value of ^ 4 x 10 ~^ for the two scales. No significant 
deviations were found related to the skewness of the 
wavelet coefficients. The analysis was repeated in the 
two Galactic hemispheres separately. This was moti- 
vated by previous findings of asymmetries related to the 
genus [QT and the N-point correlation function [42 . It 
was found that whereas there was not deviation on the 
northern hemisphere, the excess of kurtosis was even 
more remarkable in the southern region. In partic- 
ular, at the scale of 250 arcmin, the kurtosis of the 
SMHW coefficients was associated with a p-value of 
^ 2x 10~^. Again, no deviation on the skewness was no- 
ticed. Hence, these analyses indicated that the source 
for the deviation on the kurtosis of the wavelet coeffi- 
cient was related to feature/s with a typical length of 
around 10° in the sky, and located in the southern hemi- 
sphere. Thanks to the frequency /spatial properties of 
the SMHW decomposition, it was possible to study in 
detail such features. In the left panel of figure |3) the 
cleaned CMB map obtained from the fourth release of 
the WMAP data is given. In particular, I plot the opti- 
mal combination for cosmological analysis (hereinafter, 
the QVW map), obtained as a noise- weighted linear 
combination of the Q-, V-, and W-bands, previously 
cleaned via a template fitting (e.g. [50]). Regions highly 
contaminated by Galactic foregrounds, nearby clusters, 
and extragalactic point sources, have been masked. In 
the right panel, the SMHW coefficients of the previous 
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Figure 4: Right panel: cold area of the SMHW coefficients (at R = 300 arcmin), as a function of the threshold (zA 
Left panel: cold area (A^^) of the SMHW coefficients (at u = San)^ as a function of the scale (R). As in figure 
the red, green, and magenta regions represent the 32%, 5% and 1% acceptance intervals, respectively. These plots 
correspond to the analysis done by [28] . 



map are represented, at the scale R = 250 arcmin. It 
is evident the presence of the Cold Spot in the south- 
west side of the image. In particular, the centre of the 
Cold Spot is estimated to be 6> = 147° and = 209°. 
The study of the Cold Spot, through the application 
of follow-up tests, provided further evidences for its 
anomalous nature. Let me review these tests in the 
following subsections. 

2.2.2 The amplitude 

One of the most trivial statistics to study extreme val- 
ues (as the Cold Spot) in a random sample is the 
largest /smallest observation. In [118 it was established 
that the temperature of the Cold Spot was -4.57 times 
the dispersion of the SMHW coefficients at R = 250 ar- 
cmin. This cold value represented a p- value of 0.01 
(relative to Monte Carlo simulations). A more ro- 
bust statistic related to the extreme values is the MAX 
statistic, understood as the largest observation (in ab- 
solute value). For the particular case of the SMHW 
coefficients, MAX is defined, at scale i?, as: 



MAXR = max{|^(^,,(/),;i?)|}. 



(8) 



The MAX statistic is more robust than selecting the 
coldest of the extrema, since the selection of the lowest 
values could be seen as an a posteriori selection. This 
statistic was studied in [30 , showing that the Cold Spot 
was always the maximum absolute observation of the 
WMAP data at scales around 300 arcmin, representing 
an upper tail probability of 0.38% (relative to Monte 
Carlo simulations) . This value was less significant than 



the one mentioned in the previous subsection. The rea- 
son for this change is, as commented, that the MAX 
statistic is more robust than simply selecting the small- 
est values of the observations. 

2.2.3 The area 

The area above or below a given threshold is one of the 
most common statistics used to characterize the prop- 
erties of a random field. In particular, the area is the 
most commonly Minkowski functional used in the lit- 
eratur^(see, for instance, [23 ) l53l [TUB] ) . Generalizing 
this concept to the case of the wavelet coefficients, we 
can define cold {A^^) and hot (A^^) areas, at a given 
threshold u and a given scale R^ as: 



^7 = #{w{ei,4>uR)<-iy}, 

= #{wi0i,^i;R)>+,y}, 



(9) 



where the number operator #{condition^} indicates 
how many times condition^ is satisfied, for i rang- 
ing from 1 to Npix(i^). The cold and hot areas of the 
WMAP data were analyzed by [28 . It was reported 
that, whereas the hot area was consistent with the ex- 
pected behaviour for the standard Gaussian model (at 
all the scales R and thresholds u) , the cold area was not 
compatible. In particular, deviations from Gaussian- 
ity were found, again, at SMHW scales of R ^ 300°. 
The deviation took place for thresholds equal or smaller 
than —3aR (see figure |4|. The analysis per different re- 
gions of the sky confirmed that the anomaly on the 
cold area was localized in the southern-west Galactic 



-•^For 2D images, there are three Minkowski functionals, namely the contour or length, the area, and the genus. These three quantities 
are defined above/below a given threshold. 
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quadrant of the sky, and that the Cold Spot was re- 
sponsible for this anomaly. In particular, the cold area 
of the WMAP data (at the mentioned scales, and be- 
low a threshold of — 3<Ji^) was found anomalous with a 
probability of ^ 99.7%, whereas it became fully com- 
patible with Gaussian simulations, once the Cold Spot 
was not considered in the analysis. As for the case of 
the MAX statistic, a more conservative estimator (i.e., 
less dependent from the fact that the Cold Spot is nega- 
tive) can be considered, just by selecting the maximum 
value of the previous cold and hot areas: 



(10) 



This new statistic was used by [30 , finding, again, 
that the WMAP data was anomalous about thresholds 
larger than 1 3(7^^1, for scales of the SMHW of around 
300 arcmin. 

2.2.4 The Higher criticism 

Higher criticism (HC) is a relatively new statistic intro- 
duced in 2004 by [39] , and firstly applied to the context 
of probing the Gaussianity of the CMB only a year after 
by [20 . Although there is not a unique definition for the 
HC, all the forms proposed in the literature satisfy the 
same key concept: HC is a measurement of the distance 
between a given sample of n elements to a Gaussian 
probability density distribution, established by means 
of the difference between the p- value pi of a given obser- 
vation Xi — assuming it comes from a N(0, 1) — , and its 
cardinal position on the sorted list (in increasing order) 
of p-values Pi (i.e., Pi^i < Pi < Pi+i, V i = 1, . . . ,n). 
The HC associated with the sample is just defined as 
the largest value of such differences. 

This concept can be applied to the SMHW coeffi- 
cients of a given signal (e.g., the QVW map) at a given 
scale R. This was the analysis proposed by [20 . Let 
us adopt the following definition for the HC associated 
with Npix(i?) wavelet coefficients w {Oi^(j)i] R)^ at scale 
R: 

HCn,,.(R) = maxji/C^^^Ji?)} , (11) 

where the maximization is made over the quantity 
HCl^ .^(i?), that provides the difference between the 
experimental probability of the wavelet coefficients 
w {Oi^ R) at scale R and the corresponding theoreti- 
cal p-value. Such quantity reads as: 



Npix(fl) 



\/p.{R){i-Pi{R)y 



(12) 



where the p-value is given by Pi{R) = P{|N(0, 1)| > 
\w {6i^(j)i; R) \ }. The " operator indicates that the 
Npix(i^) SMHW coefficients at the scale R have been 
transformed into a zero mean and unit variance sam- 
ple. Let me remark that the p-values have been sorted 
in increasing order. 
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Figure 5: The HC^ . values obtained from the anal- 
ysis of the QVW map. The solid red line corresponds 



to the application of the equation 12 to the CMB map 
in the real space (i.e., = 0), whereas the dot-dash 
blue line corresponds to the analysis of the SMHW co- 
efficients at scale R = 300 arcmin. HCL curves are 

-L^ pix 

normalized by their respective minimum values. 



As explained in [20^, the HC represents, under cer- 
tain conditions, some advantages with respect to more 
traditional statistics designed to study the Gaussian- 
ity of a given sample. In particular, HC seems to be a 
better estimator than the MAX statistic: whereas the 
latter is designed to capture Gaussian deviations caused 
by very large values of the distribution, the HC is also 
sensitive to anomalies produced by moderate values. In 
addition, the HC can identify which values (in a given 
sample) are the ones that differ from the theoretical 
Gaussian distribution. 

In [20 it was reported that the HCn was 
above the 1% acceptance interval, again, at the SMHW 
scale R = 300 arcmin. They found that, actually, all 
the SMHW coefficients associated with HC^^,^{R) val- 
ues above the 1% acceptance interval set by CMB Gaus- 
sian simulations, where localized in the position of the 
Cold Spot. This extra test was an additional support 
to the anomalous nature of the WMAP data, and of the 
Cold Spot in particular. Results were confirmed by [3^ 
for the analysis of the second WMAP data release, re- 
porting an upper tail probability even lower than for 
the 1-year data. 

As an illustration of the HC statistic, in figure [5] 
I represent (in solid red) the HC^^^,^ values obtained 
from the QVW map in the real space, and (in dot- 
dash blue) the corresponding curve for R = 300 arcmin. 
These quantities are normalized to their minimum val- 
ues, for a better comparison. They are represented in 
an increasing order. It is remarkable that, for the case 
of the analyses performed on the SMHW coefficients, 
there is a tail of very large values of HC^ . , that are 
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Figure 6: Maps of HC^ . obtained from the analysis of the QVW map derived from the WMAP data. Left panel 
corresponds to the study of the real space case, while right panel shows the outcome of the analysis of the SMHW 
coefficients at R = 300 arcmin. Whereas for the former there are not particular signatures in the map, the wavelet 
analysis shows some prominent features, being the Cold Spot the most pronounced one. 



not present for the real space case. This is due to the 
ability of the SMHW transform to enhance features of 
a given scale and shape. In figure |6| these values are 
represented on the celestial sphere (left panel for the 
real space case, and right panel for the SMHW coeffi- 
cients). The figure indicates that there are not partic- 
ular signatures in the real space, whereas the SMHW 
coefficients at = 300 arcmin allow us for a clear iden- 
tification of the features causing the anomalous values 
of the HC^^.^{R). In particular, the key role played by 
the Cold Spot is highlighted. 

3 The characteristics of the Cold Spot 

In this Section, I summarize briefly some of the most 
important properties of the Cold Spot. I will focus in 
two major aspects: its morphology and its frequency 
dependence. 

The morphological properties of the Cold Spot are 
different depending whether we do the analysis in the 
real or wavelet space. As it was pointed out in 
the region associated with the Cold Spot, in the real 
space, appears as formed by several small cold spots. 
The amplitude of the most prominent of these spots is 
< —350/j.K with a size of ~ 1°. None of these struc- 
tures is particularly anomalous. An image of the Cold 
Spot in the real space is shown in the left panel of fig- 
ure j?! It is, however, in the wavelet space where the 
Cold Spot appears more interesting. In the right panel 
of figure [7| I present a close view of the Cold Spot af- 
ter convolution with the SMHW at a scale R = 250 
arcmin. Besides all the anomalous characteristics pre- 
viously discussed (i.e., area, HC, and MAX), the Cold 
Spot appears as a very symmetric feature. However, 
this effect could be biased since, after all, the SMHW is 



an isotropic filter and, therefore, the symmetric features 
of the Cold Spot could be amplified, erasing any possi- 
ble intrinsic anisotropy. This issue was studied in detail 
by [29 . Instead of applying an isotropic wavelet, the 
anisotropic Elliptical Spherical Mexican Hat Wavelet 
(ESMHW) was adopted. The cleaned CMB map de- 
rived form the WMAP was transformed into ESMHW 
coefficients (at the scales for which the WMAP data 
appeared as anomalous), for different ratios ( between 
the smallest and the largest axes of the ESMHW, and 
for different orientations. This work proved that the 
maximum matching between the Cold Spot and the 
ESMHW took place when ( e [0.875, 1] and, hence, 
indicated that the Cold Spot structure was quite close 
to be isotropic (assuming that the ratio of the ESMHW 
axes mimics, somehow, the ratio between the Cold Spot 
axes). 

The frequency dependence of the Cold Spot has 
been a matter of study soon after its discovery. Already 
in [118J the SMHW coefficients of the cleaned WMAP 
frequency channels (namely, Q, V and W bands) were 
computed, and the mean value of the coefficients as- 
sociated with the Cold Spot at the scale R = 250 
arcmin was estimated. No obvious frequency depen- 
dence of this mean value of the wavelet coefficients was 
found, hence, being fully consistent with the expected 
behaviour for the CMB emissiorj^ and, therefore, quite 
different to the typical frequency dependence of the 
Galactic foregrounds. 

4 The significance of the detection: the 
a posteriori issue 

One of the most questioned aspects of the WMAP 
anomalies in general, and of the Cold Spot in partic- 

^The analyzed WMAP frequency channels were in thermodynamical temperature and, therefore, the CMB appears as a frequency 
independent emitter. Notice that, since the SMHW transform is a linear operation, the same behaviour is expected for the wavelet 
coefficients. 




Figure 7: Image of the Cold Spot {0 = — 147°,^ = 209°) in the real (left panel) and wavelet (right panel) spaces. 
The homogenously filled circles correspond to positions where known extragalactic point sources have been masked. 



ular, is the issue of the actual significance of the de- 
tection. This is a very important point that is inti- 
mately linked to the blind nature of all the Gaussian- 
ity/isotropy tests that led to the report of such anoma- 
lies. 

Let me review where this problem comes from: if 
many tests are performed in a given data set, it is not 
strange that some of them report some deviation from 
the null hypothesis. It is quite usual to face the follow- 
ing situation: a set of blind tests (i.e., tests that just 
challenge the compatibility of the data with a given null 
hypothesis, Hq, and not confronting such hypothesis 
with an alternative one. Hi) claim a given incompatibil- 
ity of the WMAP data. A subsequent test is performed, 
taking into account the previous finding and, usually, 
in such a way that the initial reported deviation is now 
found at higher significance. In this procedure there 
are two weak points: the first one, already mentioned, 
is to assess the probability of finding a deviation as the 
one claimed during the first step, taking into account 
all the possible tests that were performed. The second 
one is the credibility of the probability for the follow-up 
test, where a particularity was studied in greater detail. 

As mentioned above, this is a common situation for 
the WMAP anomalies works and, therefore, the Cold 
Spot is not an exception. Several tests were made in the 
first work by [TT8 , namely, the estimation of the skew- 
ness and the kurtosis at several scales of the SMHW. 
A particular deviation was highlighted: the excess of 
kurtosis at several scales around R = 250 arcmin. Af- 
ter that, the Cold Spot was identified as a prominent 
feature, and further tests (the MAX, the cold area, the 



HC) were applied. I believe that there is not a unique 
and clear way to solve these ambiguities and, to my 
view, this point is not usually addressed in the litera- 
ture. However, whereas for the latter aspect (i.e., the 
significance for the follow-up tests) the solution is hard, 
for the former there could be some possible getaway, 
at least depending on the complexity of the prelimi- 
nary analysis. Actually, the Cold Spot is one of the 
few WMAP anomalies where this particular aspect has 
been considered with deeper interest. In fact, it was 
the matter of several papers [791 US] and, in particular, 
of [30^. 

In this last work, the significance of the first detec- 
tion was addressed, focusing in the a posteriori selec- 
tion of the statistic (the kurtosis) and the scale range 
(^ 250 arcmin). A conservative procedure to establish 
the p-value of the non-Gaussian detection, based on 
the characteristics of the analysis, was proposed. More 
specifically, since 30 statistics were applied to the QVW 
map (i.e., the skewness and the kurtosis of the SMHW 
coefficients at 15 scales), and only 3 out of these 30 
statistics were found as anomalous (i.e., the kurtosis of 
the SMHW at scales R = 200,250, and, 300 arcmin 
were outside the 1% acceptance interval — see [30] for 
details), then, it was decided to estimate the signifi- 
cance of the non-Gaussianity detection by exploring in 
how many out of 10,000 CMB Gaussian simulations it 
was observed that the skewness or the kurtosis of the 
SMHW coefficients were outside the 1% acceptance in- 
terval, at least, at three scales. The p-value obtained in 
this manner was 0.0185. This p-value can serve, as ex- 
plained in [30^, as a conservative probability related to 
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the non-Gaussianity associated with the SMHW anal- 
ysis. 

In this spirit, the fohowing up tests (e.g., the am- 
phtude, the area or the HC) can be just seen as addi- 
tional probes to explore/understand the previous devi- 
ation, rather than as independent sources for establish- 
ing a proper significance level for the detection. 

Recently, [130^ has questioned the non-Gaussianity 
found by |118j , since the excess of the kurtosis was 
clearly found with the SMHW, but it was not the case 
with other analyzing kernels (proposed in [130]), as the 
top-hat and the Gaussian filters. The authors argued 
that these tools are more natural than a wavelet like the 
SMHW and that, therefore, the selection of the SMHW 
is somehow a posteriori. Contrary to this reasoning, I 
found that the results obtained by [130 imply a dif- 
ferent conclusion: the lack of detection when analyzing 
with the top-hat and the Gaussian filter is a proof of 
the issue discussed in Section [2?T| namely, that any fil- 
tering kernel is not necessary suitable for the detection 
of any non-Gaussian feature. It is clear that some fea- 
tures (like point sources, cosmic strings or textures) are 
much better detected after applying optimal or targeted 
filters, rather than general ones (like the top-hat or the 
Gaussian functions). 

The reason why a compensate filter as the SMHW 
gets a much larger amplification as compared to uncom- 
pensated kernels as the previous ones is that it is much 
more efficient to remove the background ffuctuations 
above and below a given scale interval. Even more, it 
can be shown (e.g., [6l]) that the SMHW is close to 
the optimal or Matched filter to detect objects with 
a Gaussian-like profile embedded in CMB-dominated 
background (as it is the case in the region of the Cold 
Spot). The form of a Matched filter designed for a given 
situation is defined, not only by the shape of the fea- 
ture to be detected, but also by the statistical proper- 
ties of the background. In particular, the Matched filter 
(defined in the Harmonic space) is proportional to the 
shape of the feature, and inversely proportional to the 
angular power spectrum of the background. Therefore, 
bearing in mind that (at degree scales) the CMB is 
well described by an angular power spectrum close to 
Ci (X it is trivial to show that the SMHW is near 
to be an optimal tool for detecting a feature described 
by a Gaussian-like profile, and embedded in such back- 
ground. 

Let me remark that, as it was said in Section [2]T] 
the selection of wavelets as a suitable tool for non- 
Gaussian analysis, as it was shown by several authors in 
the past (e.g., [62] |2l [78] ) , can not be considered as an 
a posteriori choice, but, rather, as a natural option for 
studying scale-dependent phomena. In particular, the 
compensation property satisfied by wavelets make them 
extremely good analyzing kernels to amplify certain fea- 
tures, precisely because it assures a strong suppression 



of the large scale fiuctuations of the background. 
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Figure 8: Frequency dependence of the excess of the 
kurtosis of the SMHW coefficients. The plot shows the 
variation of the kurtosis as a function of the scale Kji for 
cleaned CMB maps obtained from the WMAP data at 
three cosmological frequencies (Q-,V-, and W-bands), 
and for the optimal CMB cleaned map provided by 
WMAP, as noise- weighted combinations of the previous 
maps (the QVW map). The coloured regions represent, 
as for other ffgures, the 32%, 5% and 1% acceptance in- 
tervals provided by simulations. 

Finally, it is worth recalling that the Cold 
Spot has been identiffed as an anomalous feature 
by other tools different from the SMHW: by direc- 
tional wavelets [79], scalar indices [98l [99], steer- 
able wavelets |120j , needlets [SJ, and the Kolmogorov 
stochastic parameter [56] . 

5 Some possible sources to explain the 
Cold Spot 

To find an explanation for the non-Gaussianity devi- 
ation associated with the Cold Spot is the next step, 
once its anomalous nature (i.e., non-compatible with 
the standard infiationary scenario) is accepted. With 
this aim, many efforts have been done in the last years, 
considering different sources for the observed anomaly. 
The possible causes addressed so far account for system- 
atics effects, mostly due to instrumental aspects that 
are not well understood/modelled; spurious emissions 
due to foregrounds or contaminants of the cosmological 
signal; non-accounted secondary anisotropics induced 
on the CMB photons, as the interaction with the ion- 
ized medium (e.g., the Sunyaev-Zeldovich effect) or the 
non-linear evolution of the gravitational potential (e.g., 
the Rees-Sciama effect); and, of course alternative (or 
complementary) models to the standard inffationary 
scenario (as cosmic defects). In the following subsec- 
tion, and also in Section [6) I address these possibilities. 
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Figure 9: Null tests performed on difference maps, free of any foreground and cosmological signal. The left panel 
represents the kurtosis of the SMHW coefficients, as a function of the scale i?, for the Ql — Q2 map, i.e., for a map 
built as the difference of the difference assemblies observations at 41 GHz. The middle panel is as the previous one, 
but for the VI — V2 difference map, i.e., at 61 GHz. Finally, the right panel provides the variation of the SMHW 
coefficients as a function of the wavelet scale, but for the difference map obtained with the 4 difference assemblies 
of WMAP at 90 GHz, i.e., Wl — W2 + W3 — W4. The blue asterisks represent values for the data, whereas red, 
green and magenta regions represented the acceptance intervals at 32%, 5% and 1%, respectively, determined by 
simulations. 



starting from those hypotheses that are less dramatic 
from the point of view of strong implications for the 
standard cosmological model. 

5.1 Systematics 

To test the influence of unknown/unmodeled system- 
atics on the non-Gaussianity deviation is, as one could 
imagine, a very hard task. The analyses that can be 
carried out to probe such sources are, basically of two 
types. 

One of these analyses are consistency tests. As 
most of the CMB experiments, WMAP satellite can 
provided us with cleaned CMB maps (for instance, fol- 
lowing the template fitting approach described by [50] ) 
for several detectors. Therefore, focusing in the non- 
Gaussian deviation associated with the Cold Spot, an 
obvious procedure would be to check whether the ap- 
plication of the different statistical tools reveals that 
such a feature is associated with only one detector, or 
a smaller set of detectors. If this were the case that 
would be a clear indication for a lack of consistency 
and, therefore, that the non-Gaussianity detection is 
associated with a given instrumental feature. This was 
done by |118j and [28 for the kurtosis (Kr) and the 
area {A^) of the SMHW wavelet coefficients. No incon- 
sistency was found: the excess of kurtosis and of area 
was found to be the same for every difference assem- 
bly. As an illustration, in [8] the kurtosis of the SMHW 
coefficients (as a function of the scale Kr) for 4 dif- 
ferent CMB maps is presented. In particular, results 
for the Q-, V-, and W-band cleaned CMB maps are 
shown, together with the curve obtained from the anal- 
ysis of the optimal QVW-map. These curves are quite 
similar, which indicates that the non-Gaussian signal 



is presented in all the WMAP detectors, at a similar 
level. 

The second type of analyses are null tests. The 
cosmological frequencies of the WMAP satellite (i.e., Q 
band at 41 GHz, V-band at 61 GHz and W-band at 
94 GHz) are made from more than one difference as- 
sembly. Hence, just subtracting difference assemblies 
at the corresponding band can produce noise maps per 
frequency. Neglecting small differences from the op- 
tical beams and the band-pass widths, the CMB and 
foreground emissions have been cancelled out in this 
new map. Therefore, the application of the statistical 
tools to these difference maps helps to check whether 
the non- Gaussian signal is a noisy artifact (if such sig- 
nal is still present) or not (if consistency with Monte 
Carlo simulations is found). These was done by |118j 
and [28 , again, for the kurtosis and the area of the 
SMHW coefficients, respectively. As for the previous 
type of systematics probe, there was a clear indication 
that the non- Gaussian signal was not related to any 
instrumental signature. As an example, in figure [9j 
the result obtained by [118 for Kr is presented. On 
the left panel, I represent the variation of the SMHW 
wavelet coefficients for the difference map constructed 
at 41 GHz as Ql — Q2. Similarly, the results for the 
VI — V2 map at 61 GHz are provided in the middle 
panel. Finally, in the right panel, I give the output for 
the difference map obtained at the 94 GHz band as the 
combination Wl - W2 + W3 - W4. Results for the 
WMAP data is given as blue asterisks, whereas, as for 
previous figures, the red, green, and magenta regions 
provide the 32%, 5% and 1% acceptance intervals, re- 
spectively. 

Summarizing, consitency and null tests do not 
reveal the presence of systematics behind the non- 
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Gaussianity associated with the Cold Spot . Besides 
these test, it is important to remark that the angular 
size associated with this feature is ~ 10° in the sky. It 
is not trivial to think in a systematic effect affecting at 
this scale, and providing such localized feature in the 
sky as the Cold Spot . 

Finally, the subsequent releases of the WMAP 
data (where the modelling of the instrumental prop- 
erties have been improving with time) have shown that 
there are not changes in the no n- Gaussianity deviation, 
except for a slight increasing on its significance, which 
reflects the higher signal-to-noise that WMAP data is 
getting as observational time increases. 
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5.2 Foregrounds 

Astrophysical contaminants or foregrounds are the next 
possible origin for the non- Gaussianity associated with 
the Cold Spot,. It is well known that foregrounds are 
highly non- Gaussian signals. It is worth commenting 
that, although the Cold Spot is negative, it is still pos- 
sible to think in an additive source (as foregrounds are) 
as a feasible explanation. To understand this point, it is 
important to recall that the QVW map (that, as I said 
before, is commonly adopted in the literature for cosmo- 
logical analyses) is obtained as a noise-weighted linear 
combination of cleaned CMB maps at different frequen- 
cies. These maps (at 41, 61 and 94 GHz) are produced, 
as mentioned previously, via a template fitting (see, for 
instance, [50 ). Therefore, any over-subtraction of fore- 
grounds templates could caused a foreground residual 
in the form of a cold emissiorj^ 

Since, as mentioned before, the size of the Cold 
Spot is of several degrees, it is really hard to believe 
that this feature could be associated with residuals from 
point sources (i.e., from radio and infrared galaxies). 
Therefore, only galactic emissions (as synchrotron, free- 
free, and thermal and spinning dust) could be responsi- 
ble for a large feature as the Cold Spot. However, notice 
that the Cold Spot is located at 57° from the Galactic 
plane, and at a longitude of 209°. In other words, the 
Cold Spot is placed in a region of low Galactic contam- 
ination. According to the previous reasoning, it is al- 
ready hard to make compatible the presence of the Cold 
Spot with a given Galactic emission. Even though, of 
course, the issue has been a matter of discussion. I 
review here some of these analyses. 

The most obvious test is, of course, to check 
whether there is or not any frequency dependence of 
the statistical estimators that indicated the Gaussian 
deviation. For instance, in |118( [29] , the kurtosis of 
the SMHW coefficients {Kr) was studied for the dif- 
ferent CMB cleaned maps obtained at the Q-, V-, and 
W-bands. 

^In some works, a VW map is adopted, i.e., a map built as a 



Figure 10: Frequency dependence of area of the SMHW 
coefficients. The plot shows the area of the wavelet coef- 
ficients at R = 300 arcmin for different thresholds (from 
top to bottom: 3, 3.5, and 4(7^^), as a function of the 
frequency. The circles represent the values obtained for 
cleaned CMB maps obtained from the WMAP data at 
three cosmological frequencies (Q-,V-, and W-bands). 
The solid lines represent the areas, at different thresh- 
olds, for the noise- weighted QVW map. 

The results obtained for this study are given in 
figure |8] The kurtosis Kr is presented for the 41, 61, 
and 94 GHz channels, and for the noise- weighted lineal 
combination (the QVW map). It is remarkable the high 
similarity of the curves. The pattern of the kurtosis, as 
a function of the SMHW scale, is the same for all the 
maps. The same is observed for its normalization. An 
equivalent test can be done for the area of the SMHW 
coefficients (A^), as proposed by [28l[29). Results, at 
R = 300 arcmin, are presented in figure [lO] Notice 
that the agreement of the area of the SMHW coeffi- 
cients (above threshold u =3, 3.5, 4<Ji^) among differ- 
ent frequency bands (Q, V, and W) and the combined 
QVW map (represented by the solid lines in the figure) 
is very high. These kind of tests show that there is not 
any evident frequency dependence of the statistics asso- 
ciated with the non-Gaussian deviation and, therefore, 
that such anomaly is fully consistent with the expected 
behaviour for a CMB feature. 

Additional tests supporting this idea have been 
proposed in the literature. First [118 , the kurtosis of 
the SMHW coefficients was analyzed for a CMB-free 
map, constructed as the combination of the 4 differ- 
ence assemblies at the W-band minus the sum of the 4 
ones at Q- and V-band (i.e., Wl + W2 + W3 + W4 - 
Ql - Q2 - VI - V2). This kind of map could have 
a contribution of the Galactic contaminants (outside a 
given observing mask), since foreground emissions are 
not frequency independent. This analysis did not show 
any significant deviation from the expected behaviour 

^-weighted linear combination of the cleaned V- and W-bands. 



11 



from Gaussian simulations, and, therefore, discarding 
the presence of significant foreground residuals. 

Second, as suggested by [29 , different CMB re- 
coveries from the WMAP data (where independent and 
alternative component separation approaches were fol- 
lowed) could be analyzed. In particular, the cleaned 
CMB maps obtained by [112J were probed. The kurto- 
sis, the area, and the MAX of the SMHW coefficients 
do not changed significantly from the different CMB 
maps. 

[29 explored a more complicate scenario: a sit- 
uation in which combinations of different foreground 
emissions could mimic, in the region of the Cold Spot, 
the behaviour associated with the CMB, i.e., a fre- 
quency independent global emission. To check this 
possibility, several templates were used as tracers of 
the Galactic foregrounds, namely, the Rodhes/HartRaO 
2326 MHz ^8] radio survey for synchrotron, the Hq, 
by ^48^ for the free-free, and the thermal dust model 
by [IZl- Authors studied the expected contribution of 
foregrounds in the region of the Cold Spot, and they 
found that, taking into account the uncertainties in the 
extrapolation of these templates from their original ob- 
servations to the WMAP frequency range, it was pos- 
sible to find a global Galactic emission that was nearly 
frequency- independent from 41 to 94 GHz (i.e., from 
Q- to W-bands). However, it was found that the emis- 
sion was at a level of one order of magnitude below 
the Cold Spot temperature. It was checked that, even 
accounting twice for that hypothetical foreground emis- 
sion, it was not possible to reconcile observations with 
the Gaussian model. 

All these tests on the impact of the foregrounds 
indicated that the non- Gaussian signal associated with 
the Cold Spot was fully consistent with a CMB like 
frequency dependence, and that the role played by as- 
trophysical contaminants was negligible. 



5.3 The Sunyaev-Zeldovich effect 

After checking that the possible impact of systemat- 
ics and foregrounds on the non-Gaussian detection is 
very low, the next step is to study whether secondary 
anisotropics of the CMB could explain the anomalous 
nature of the Cold Spot. 

The Sunyaev-Zeldovich effecl|^(SZ) could be a po- 
tential candidate to explain the anomaly. Two major 
reasons support this possibility: first, the size of the 
Cold Spot is nearly compatible with the fiuctuations 
caused by the nearest clusters of galaxies, and, sec- 
ond, this fiuctuations (in the frequency range covered 
by WMAP) produce cold spots [11 Oj . 
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Figure 11: Best-fit of the WMAP data —across the 8 
difference assemblies for the Q-, V-, and W- bands — to 
a CMB (solid line) and a SZ spectra (dashed line). The 
fit is performed to the minimum amplitude of the Cold 
Spot in wavelet space, at a scale of R = 300 arcmin. 
Error bars are obtained from Gaussian simulations. 

There is not evidence for the existence of any 
large cluster in the direction of the Cold Spot. How- 
ever, [15 reported the presence of a large concentra- 
tion of galaxies in that direction — referred to as the 
Eridanus super- group — , that could account for a total 
mass of ^ lO^^M© (see [32^ for details). 

This scenario implies, therefore, that the SZ could 
explain the nature of the Cold Spot, at least, partially. 
In other words, it could be possible that a combina- 
tion of a SZ contribution plus a large (but not anoma- 
lous) CMB fluctuation could account for the Cold Spot 
emission. This was studied by [29 , and the results 
are summarized in figure [TT] In this figure, the val- 
ues of the SMHW coefficients w{9^(j)]R)^ at the po- 
sition of the Cold Spot and at = 300 arcmin, are 
given for the eight CMB cleaned maps at the Q-, V-, 
and W-bands. The corresponding error bars were com- 
puted from CMB plus noise simulations, corresponding 
to the instrumental properties of these detectors, and 
convolved with the appropriate SMHW kernel. 

Three different fits to the data were explored. 
First, a pure CMB spectrum was using, giving a very 
good fit with a reduced = l-OO- Such fit is repre- 



sented by the solid line in figure [IT] Second, a pure SZ 
spectrum was tested, obtaining a very poor fit (dashed 
line) with a reduced = 9.12. Finally, a joint fit to a 
CMB plus SZ spectra was explored, obtaining an am- 
plitude for the SZ spectrum consistent with zero, and a 
reduced quite similar to the first case. These results 
confirm that the frequency dependence of the Cold Spot 
is consistent with a CMB-like spectrum. They also rule 

"^It is produced by the inverse Compton interaction of the CMB photons, as they cross the hot electron gas that is found in clusters 
of galaxies. 
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out the possibility that the SZ is playing any significant 
role. 

5.4 The late evolution of the large scale struc- 
ture 

Another secondary anisotropy that could explain the 
anomalous nature of the Cold Spot, is the one due to 
the non-linear evolution of the gravitational field: the 
so-called Rees-Sciama effect (RS). In particular, it is 
known [76l [77] that voids in the large scale structure 
could induce a negative non-linear anisotropy in the 
CMB photons. The size of such secondary anisotropics 
depends on the proper size of the void and its redshift. 
Therefore, as for the SZ, the RS is another potential 
candidate to explain the anomalous nature of the Cold 
Spot. 

Extra support for this hypothesis came from two 
different paths. On the one hand, theoretical works 
as [nUinSl proposed that a very large void (of ~ 300h~^ 
Mpc) and located at low redshift <C 1) could pro- 
duce large negative CMB fluctuations such as the Cold 
Spot, even with modest density contrast values (i.e., 
in a quasi-linear regime). On the other hand, it was 
suggested [102^ that the NVSS catalogue [24 seems to 
show, at the position of the Cold Spot, a lack in the 
number count of radio galaxies. 

Against these ideas, some criticisms can be made. 
First, this kind of voids are not observed and, even 
more, according to current N-body simulations j22] 
they are extremely rare events (~ 13cr, i.e., much more 
rare than the Cold Spot deviation itself, that was a 
98.15% event!) in the standard cosmological frame- 
work. Second, the claim made on the NVSS data has 
been recently questioned by |lQ8j , suggesting that such 
finding was an artefact caused by possible systematics 
related with the NVSS data processing, the statistical 
procedure and the a posteriori selection of the Cold 
Spot position. 

In addition, observational campaigns on the region 
of the Cold Spot were made by [54 with the MegaCam 
on the Canada-France-Hawaii Telescope, and by [12] 
with the VIMOS spectrograph on the Very Large Tele- 
scope. Both works have reported that the large scale 
structure in that direction, and up to z ~ 1, is fully 
consistent with the standard model, and that there is 
no evidence of such large voids as those required by 
the non-linear evolution of the gravitational potential 
hypothesis. 

6 A plausible explanation: cosmic tex- 
tures 

In the previous Section, I have presented an overview 
of the several works carried out to establish whether 
the non-Gaussianity detection associated with the Cold 



Spot could be explained in terms of systematics, fore- 
grounds, and secondary anisotropics as the Sunyaev- 
Zeldovich and the Rees-Sciama effects. None of these 
possibilities seem to provide a satisfactory explanation 
and, therefore, other sources should be investigated. 

In this Section, I pay attention to the suggestion 
made by [31]: the Cold Spot could be caused by a 
cosmic texture. Cosmic textures jll4j are a type of 
cosmic defects. They are supposed to be generated at 
some stage of the early Universe, associated with the 
symmetry-braking phase transitions that are predicted 
by certain theoretical models of high energy physics 
(see, for instance, [125' and references therein for a 
much more detailed explanation). In short, defects can 
be understood as space regions of a given phase state, 
surrounded by a space already in a new phase. In some 
cases, as for textures, these regions could collapse and, 
therefore, left an imprint on the CMB photons. 

Among the different types of cosmic defects, tex- 
tures are the most plausible candidate to explain the 
anomalous nature of the Cold Spot, since the interac- 
tion of the CMB photons with the time variation of 
the gravitational potential, associated with an even- 
tual collapse of the texture, produces spots in the CMB 
fiuctuations [115 . Even more, cosmic textures are ex- 
pected to be a possible source of kurtosis deviation, 
whereas the expected level of skewness is almost negli- 
gible (at least for values of the symmetry-breaking en- 
ergy scale compatible with current observations). This 
is caused by the even probability of textures producing 
cold and hot spots and, therefore, providing a nearly 
symmetric distribution of temperature fiuctuations. In 
fact, the equilateral /nl expected from textures goes as 
/nl ^ 1-5 X 10~^^^iIjq |lQ3j . where V^o is the symmetry- 
breaking energy scale, measured in Gev. For typical 
limits in ipo imposed by the CMB angular power spec- 
trum analysis (e.g., [9) I116j ). the expected equilateral 
/atl is ^ 10~^, i.e., a tiny value well below the current 
constraints [^ . 

The isotropic shape of the temperature fiuctua- 
tions related to these spots can be approximated, at 
least at small angular distances as [92] : 

AT 1 

^W = ±e^ (13) 

where i9 represents the angular distance from the centre 
of the spot, and i9c is a characteristic scale parameter 
of the spot — that is related to the redshift of the spot 
and the dynamics of the Universe [31] . The amplitude e 
is proportional to the symmetry-breaking energy scale 
ipo'. s = Sir'^GipQ. It is worth remarking that, accord- 
ing to cosmic texture models, the amplitude s is the 
same for every single spot generated by the collapsing 
defects. 
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Figure 12: Radial profile (solid blue/red lines) of a 
cosmic texture, as given in equation [M] The param- 
eters defining this profile are e = 10~^ and = 5°. 
The dashed lines represent the truncated values of the 



profiles given by equation 13 (blue), and the Gaussian 
function (red). See text for details. 

As mentioned above, the profile proposed in equa- 
tion [l3] is only valid up to small relative distances (e.g., 
^ '^c)[3 In order to have a profile valid at larger 
angular distances, [31] proposed to extend the profile 



given by equation 13 from its half-maximum, following 
a Gaussian function. The extension is done by impos- 
ing continuity, both, of the profile itself and of its first 
derivative. 

Taking into account these conditions, the assumed 
profile for the CMB temperature fluctuation caused by 
a collapsing cosmic texture would be given by: 



AT 



= ± < 



r(^'+^') 



if ^ < ^* 



if ^>^*, 



(14) 



where = >/3/2i^c- 

In figure [T2j I show the radial section of the above 
profile (solid line). The red and blue parts of the 
solid line correspond to the Gaussian function and to 



the profile of equation [13) respectively. The dashed 
lines represent the truncated regions for both curves. 
The cosmic texture parameters used in this profile are 
1? = 5° (i.e., similar to the SMHW scale at which the 
Cold Spot appears as anomalous) and e = 10~^. This 
value corresponds to a symmetry-breaking energy scale 
of i/jQ = 1.13 X 10^^ GeV, that corresponds to a con- 



servative upper limit imposed by CMB measurements 
(e.g. BDli). 

Let me remark that the results obtained with a 
profile as the one given in equation [14) and that are 
reviewed in the next Section, do not depend very much 
with the specific function adopted for the extrapolation. 
Similar results are obtained, for instance, when an ex- 
ponential function or a SMHW-like kernel are used. 

6.1 The Bayesian framework 

Attending to the issues discussed in the previous sub- 
section, a cosmic texture could be a strong candidate to 
explain the Cold Spot: textures produce spots on the 
CMB temperature fluctuations, they are non-Gaussian 
signals, and, depending on their amplitude (or the 
symmetry-breaking scale), they could be compatible 
with current constraints on the role played by cosmic 
defects on the structure formation and evolution of the 
Universe. 

[31 proposed to make use of the texture profile 



of equation [14) to perform a hypothesis test to decide 
whether the WMAP data (in the position of the Cold 
Spot) is more likely to be described by a large (but not 
anomalous) CMB spot (i.e., the null Hq hypothesis) or 
by a cosmic texture of amplitude e and size i9 added to 
a random Gaussian and isotropic CMB field (i.e., the 
alternative Hi hypothesis). 

The optimal way of performing such hypotheses 
test is within the Bayesian framework. Bayes' theo- 
rem states that, given a data set D and some unknown 
parameters O (defining a given model in the context 
of a given hypothesis Hi), the posterior probability of 
the parameters/model given the data P(0 | D^Hi), is 
related to the likelihood P{D \ Q^Hi) (i.e., the proba- 
bility of the data given the parameters/model) as: 

p{D I e,Hi)P{e I Hi) 



p{e\D,Hi) = 



(15) 



P{D I H{) 

where P(9 | Hi) is a measurement of our a priori 
knowledge about the parameters/model (i.e., the prior), 
and P{D \ Hi) is a constant (i.e., it does not depend on 
the parameters/model) called Bayesian evidence (BE). 
The BE is nothing but the average likelihood with re- 
spect to the prior, 

P{D I H,) = j P{D I e, H,)P{e I H,) d6, (16) 

and it is a largely used mechanism to perform hypothe- 
ses test. In particular, its role on different cosmology 
fields has been quite remarkable during the last years 
(e.g., [86j for dark energy studies, anisotropic 
models of the Universe expansion, |87] for studying dif- 
ferent re- ionization models, [T7] for point source detec- 
tion, and [122 for exploring non-standard inflationary 
models). 



^To my knowledge, there is not any analytical or numerical solution for the full profile of a given cosmic texture yet. 



14 



The importance of BE for hypotheses test is clear. 
First, it is obvious that the quantity that we would like 
to obtain is a measurement of the probability of a given 
hypothesis i^o, given the data, i.e., P{Hq \ D). This 
probability can be written, attending to the probability 
multiplication rule^ as: 



P{Ho I D) = P{D I Ho) 



P{D)' 



(17) 



i.e., it is proportional to the BE and to the probabil- 
ity of the hypothesis, and inversely proportional to the 
probability of the data. Under certain circumstances, 
the probability of the hypothesis could be known, but, 
however, it is not the case for the probability of the 
data. In other words, we only can learn about the prob- 
ability of the hypothesis Hq up to a factor. Therefore, 
what we can extract is a relative measurement of the 
probability of two hypotheses {Hq and i^i), given the 
same data set D. This relative measurement is called 
the posterior probability ratio^ p, and reads: 



P{H^ I D) El P{Hi 



P{Ho\D) EoP{Ho)' 



(18) 



where, for simplicity, I re- write the BE, P{D \ Hi)^ as 
Ei. Hence, if p > 1, we can conclude that the hypoth- 
esis Hi is favoured by the data with respect to Hq. In 
some cases, there is not a clear choice for the probabil- 
ity of the hypotheses. In this case, empirical rules for 
the ratio of evidences — as the Jeffreys' scale [67] — are 
usually adopted. 

Therefore, the procedure required to explore 
whether the Cold Spot is more likely to be explained 
in terms of a cosmic texture {Hi) rather than by a 
Gaussian CMB fluctuation (i^o) is clear: first, the like- 
lihood is computed for both hypotheses; second, the 
BE is estimated, taking into account the adequate pri- 



ors (equation 16); finally, the posterior probability ratio 
(equation 18) is evaluated, making use of suitable a pri- 
ori probabilities for Hq and Hi. 

This was the procedure followed in [31 . It is 
straightforward to show that, since the noise term is 
caused by standard CMB Gaussian fluctuations and 
instrumental Gaussian noise, the likelihood function 
reads as: 



P(I)|e,i70^exp(-xV2), 



(19) 



where = {D - T(6)) C'^ {D - T(6))^. The corre- 
lation matrix C accounts for the full Gaussian CMB 
and noise correlations — i.e., C = S -\- where Sij oc 
C £ Pi {cos 6 ij) and riij = cfdij^ being the Legen- 
dre polynomials, 6ij the angular distance between the 
pixels i and j, cr| the instrumental noise contribution 
to pixel i, and Sij represents the Kronecker delta. D 
represents the data (i.e., the QVW map) and the func- 
tion T(B) represents the model behind the hypotheses 



— i.e., equation 14 for Hi, and = for Hq. Finally, ^ 
denotes standard matrix transpose. 

The priors adopted by [31] for the parameters 
= (e , i?c) were chosen attending to observational con- 
straints and cosmic texture simulations. In particular, 
the prior on the amplitude was < 10~^, whereas the 
prior of the size 1° < '^^c ^ 15° was assumed. The am- 
plitude prior was uniformly distributed, and, as men- 
tioned above, it is a conservative constraint imposed 
from the contribution of cosmic defects to the CMB 
angular power spectrum. The size i^c follows a scale- 
invariant law, and the limits comes form texture sim- 
ulations. Textures below 1° should be smeared out by 
photon diffusion and, in addition, they would be related 
to collapsing events above redshift ~ 1000, which would 
not affect the CMB image. The upper limit is due to the 
unlikely probability of generating such large textures in 
the finite celestial sphere. Even so, [31 tested that re- 
sults were not specially sensitivity to the prior selection, 
since the likelihood was clearly peaked, within a region 
of the parameter space clearly allowed by observations 
and texture models. 

The marginalization of the posterior probability 
in equation [Tsj led to the determination of the tex- 
ture parameters, obtaining e = 7.3^3 g x 10~^ and 

i^c = 4.9°t2;4o at 95% confidence. The BE ratio was 
150, which, in terms of the empirical rules [67 , is a 
strong indication that the texture hypothesis for the 
Cold Spot is favoured over the isotropic and Gaussian 
CMB fluctuation option. Adopting a ratio for the prob- 
ability of the hypotheses given by the fraction of the 
sky that is covered by a cosmic texture as large as the 
one required for the Cold Spot (^ 0.017), the posterior 
probability ratio was p = 2.5, which also favours the 
texture hypothesis. 

It is worth mentioning that the estimated value 
for the texture amplitude could be affected by selection 
bias. In [31 it is established that such bias could pro- 
vide an overestimate of the texture amplitude by a fac- 
tor of 2. This bias is caused because the texture ampli- 
tude is estimated in a low signal-to-noise regime, were 
the features placed in large background fluctuations are 
more easily detected. Even so, the estimated value for 
the texture amplitude, e = 7.3 x 10~^, would imply 
a symmetry-breaking energy scale of tjjo = S.7 x 10^^ 
GeV, which, on the one hand, is fully compatible with 
more recent constraints imposed from the analysis of 
the CMB temperature and polarization power spectra 
(e.g. [116j ). and, on the other hand, is in agreement 
with the predictions of most of the models for parti- 
cle physics. In addition, by relating the angular size of 
the CMB texture profile to the cosmological parameters 
defining the geometry and evolution of the Universe, it 
was possible to establish that the texture collapse (that 
generated the CMB profile in equation 14 ) occurred at 
redshift z ^ 6. 
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Figure 13: Left panel, CMB temperature fluctuations obtained from WMAP data, centred in the position of the 
Cold Spot. Middle panel, best-flt of a texture proflle to the data, according to equation 14 Right panel: CMB map, 
at the position of the Cold Spot, after the subtraction of the best-flt texture model. 



The flgure 13 shows the effect of correcting the 
QVW map from the texture emission. In the left panel 
the region of the sky where the Cold Spot was identifled 
is shown. The middle panel presents the best-flt of the 
texture proflle, according to the parameters previously 
mentioned. On the right panel, I present the resulting 
map after the subtraction of the estimated cosmic tex- 
ture contribution. The Cold Spot is noticeably reduced. 
More quantitative measurements were made by [31 . It 
was proved that, if cosmic textures (adequate to the pa- 
rameters flxed by the Cold Spot analysis) were added 
to isotropic and Gaussian simulations of CMB signals, 
as seen in the QVW WMAP map, then the kurtosis 
of the SMHW became compatible. Even more, [M] 
showed that the WMAP map (corrected from the cos- 
mic texture contribution in the location of the Cold 
Spot) was not compatible with anisotropic patterns for 
non-standard expansions of the Universe (in particular, 
for Bianchi Vllh models), as it was previously the case 
for the uncorrected data (e.g. j66l 113]). 

Finally, let me remark that [32^ proposed a sim- 
ilar approach as the one described in this Section to 
study (attending to spatial templates) whether SZ and 
RS could provided more suitable hypotheses than the 
standard isotropic and Gaussian model. This Bayesian 
analysis indicated that neither of these hypotheses is 
favoured. 



6.2 Follow-up tests 

The studies described in Section |5] and in the previous 
subsection, clearly indicates that, among all the realis- 
tic sources that could explain the anomalous nature of 
the Cold Spot, only the cosmic texture hypothesis re- 
mains as a feasible option. The results obtained by |31] , 



and reviewed in the previous subsection, have to be un- 
derstood as a (clear) indication that the cosmic texture 
is plausible. However, before accepting it as the final 
explanation, it should be conflrmed by additional tests. 
In particular, if the cosmic texture hypothesis would be 
the right one, then there are some clear predictions that 
could be tested (at least in the near feature, once on- 
going/upcoming experiments as SPT, ACT, QUIJOTE 
and ALMA are fully operative). 

In this subsection I comment on the three most 
obvious follow-up tests that could help to discard, or 
accept, the cosmic texture. These follow-up tests are: 
the searching of more textures, the local polarization of 
the CMB, and the local CMB lensing. These foreseen 
effects were flrstly pointed out by [31J. 

6.2.1 Looking for more textures 

If a texture were found in the location of the Cold Spot, 
then, attending to cosmic texture models, there should 
be more cold and hot spots randomly distributed across 
the sky. In fact, the distribution of CMB spots caused 
by cosmic textures follows a scale- invariant law: 



Nsp{> ^c) 



(20) 



i.e., the number of cold/hot spots with a scale equal or 
larger than ^c-, is inversely proportional to (see [109j 
for the grounds and [31 for a derivation). The v and n 
parameters are associated with the physics of the cos- 
mic texture models, and their speciflc meaning is out 
of the scope of this review. Let me remark here that, 
according to state-of-the-art simulations, these param- 
eters are well consistent with v ^2 and n ^ values, 
when is expressed in radians. Before continuing the 
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discussion, let me stress that the number of expected 
cold/hot spots with a scale i^c equal or larger than 5° 
(i.e., as the Cold Spot scale) is ~ 1. In other words, 
the fact that we only found a Cold Spot of a scale sim- 
ilar or larger than the Cold Spot size is fully consistent 
with the cosmic texture scenario. This fact, of course, 
is an extra support for the texture hypothesis causing 
the Cold Spot emission. 

Let me come back to equation 20 A straightfor- 
ward calculation tells as that, if the cosmic texture hy- 
pothesis is correct, then the CMB temperature fluctu- 
ations should contain ^ 28 spots with a scale i9c > 1°, 
and ~ 7 spots with a scale i9c > 2°. Poisson errors 
can be safely assumed for these numbers. Therefore, 
the number of cold/hot spots expected in the WMAP 
data are 1 < 7Vsp(> 2°) < 13 at the 95% confldence 
level. Current work is in progress to check this pre- 
diction, by using a fast cluster nesting sampling algo- 
rithm — MULTINEST [45 — to explore the posterior 
probability ratio. There are well founded hopes to flnd 
new textures. In particular, some non-Gaussian anal- 
yses as |12Q[ [94 l |56 reported some hot/cold spots (in 
addition to the Cold Spot) as potentially anomalous. 



6.2.2 The polarization of the CMB 

The effect of a collapsing texture on the passing by 
CMB photons is nothing but a secondary anisotropy 
of the CMB fluctuations, whose origin is merely grav- 
itational. Hence, the effect of such gravitational phe- 
nomenon on the E-mode polarization is almost negli- 
gible (only vector modes would be affected, which are 
well below the scalar mode contribution). 

Strictly speaking, this lack of polarization is not 
a unique signature produced by cosmic textures. As 
I said, any secondary anisotropy of gravitational ori- 
gin would cause it. However, these other effects (as 
huge voids) are quite implausible explanations (see Sec- 
tion |5.4[ ). For that reason, this effect is a valid follow-up 
test to probe the texture hypothesis. 

The procedure is simple: to compare the E-mode 
polarization in the position of a temperature spot (as 
large and extreme as the Cold Spot is), under two dif- 
ferent hypotheses, the null or Hq one (i.e., the temper- 
ature spot is caused by a Gaussian fluctuation) and the 
alternative or Hi option (i.e., the temperature spot is a 
secondary anisotropy caused by the collapse of an evolv- 
ing texture). In fact, as proposed by |124j . the best dis- 
criminating measurement is the T-E correlation, rather 
than simply the E signal. This cross-correlation is ex- 
pected to be close to zero for the Hi hypothesis. The 
approach suggested in this work was to estimate the 
correlation of the T and E proflles around the position 
of the spot temperature signal. This statistic was com- 
puted for many simulations according to the Hq and Hi 
hypotheses, and a hypothesis test was performed, via 



the deflnition of an optimal Fisher discriminant statis- 
tic (e.g. [27]). 
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Figure 14: The curve indicates the probability (the 
p- value) for rejecting the texture hypothesis, attending 
to the T-E correlation, as a function of the instrumen- 
tal noise sensitivity in polarization. From left to right, 
the vertical lines indicate the noise level associated with 
QUIJOTE, Planck, and WMAP experiments. 

The method was applied to probe the capabilities 
of current and upcoming CMB experiments for discrim- 
inating between the two hypotheses. In particular, the 
cases of WMAP, Planck [m], and QUIJOTE [100] were 
considered. 

As the major conclusion, it can be established that 
the discrimination power of the T-E correlation is not 
very high. In fact, for an ideal noise- free experiment, 
at a power of the test of 0.5, the signiflcance level is 
(up most) 0.8%. The reason for this limitation is that 
the characteristic size of the Cold Spot is ~ 10°, which, 
roughly, corresponds to multipoles of ^ ~ 40. It hap- 
pens that, at this multipole scales, the T-E angular 
cross-power spectrum is very close to zero, already for 
Hq and, therefore, it is hard to discriminate it from Hi. 

The results are graphically summarized in flg- 
ure 14 I represent the signiflcance level or p- value (for a 
power of the test of 0.5) as a function of the instrumen- 
tal noise sensitivity for the E-mode polarization. The 
three experiments previously mentioned are indicated 
as vertical lines. Notice that the WMAP instrumental 
characteristics do not allow for any signiflcant discrim- 
ination between hypotheses. With Planck, it would be 
possible to reach a modest 7% signiflcance level. Fi- 
nally, the upcoming QUIJOTE experiment would al- 
low for a 1.4% detection. However, if the signiflcance 
of the Cold Spot obtained via the flrst non-Gaussianity 
wavelet analysis (i.e., 1.85%, as mentioned in Section [4| 
is also considered, [124 claims that a joint T and T-E 
signiflcance detections of 0.025% and 0.12% could be 
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imposed on the texture hypothesis, for the QUIJOTE 
and Planck experiments respectively. 

6.2.3 Gravitational lensing 

Besides the lack of polarization discussed in the previ- 
ous subsection, the lensing of the CMB photons would 
be another foreseen effect caused by the gravitational 
field generated by the collapsing texture. This point 
was recently addressed by [35j. This work studied the 
capabilities of small-scale CMB experiments as ACT to 
detect a possible lensing effect occurring at the position 
of the Cold Spot, caused by the gravitational field of a 
texture placed at 2; = 6, and with a typical scale of 
5°, associated with a symmetry-breaking energy scale 
of ipQ = 4.5 X 10^^ GeV (i.e., comparable to the param- 
eters determined by [31]). 

As for the polarization test, the power of the lens- 
ing analysis to probe the existence of the texture is 
relatively modest. For instance, detection is made at 
the 3<j level after 1000 minutes of integration time. In 
other words, this test would require of a dedicated ob- 
servational campaign. 

7 Conclusions 

In this article I have presented a comprehensive 
overview of the Cold Spot. Since its detection in 2003 
by [118 , this feature proved to be one of the most in- 
triguing anomalies found in the WMAP data. 

The Cold Spot was detected after performing a 
non-Gaussianity test on a cleaned CMB (obtained via 
a template fitting of the WMAP difference assemblies 
at Q-, V-, and W-bands). The non-Gaussian analy- 
sis was performed by comparing the values obtained 
for the skewness and the kurtosis of the SMHW coeffi- 
cients (at several scales) for such cleaned CMB map to 
the distribution expected from isotropic and Gaussian 
CMB realizations. This analysis indicated an excess of 
the kurtosis at SMHW scales R of around 300 arcmin. 
Subsequent analysis of the SMHW coefficients, based 
on the area above/below a given threshold, the Higher 
Criticism^ and the maximum value, agreed in detect- 
ing the non-Gaussian deviation occurring at the same 
scales, and confirmed the peculiar role played by the 
Cold Spot. Avoiding any possible a posteriori choice of 
statistics, a conservative significance detection level of 
1.85% was placed by [30 . 

The Cold Spot was found to be highly 
isotropic [29^, and the impact of possible systemat- 
ics and residual foreground contamination was dis- 
carded [118[ [28| [29| [30] . Equivalently, some secondary 
anisotropics (potentially responsible for the Cold Spot 
emission) as the Sunyaev-Zeldovich and the Rees- 
Sciama effects, that, in principle, could be accepted 
as valid solutions, were show to be unlikely explana- 



tions [211 [32]. 

A plausible explanation in terms of a cosmic de- 
fect was addressed by [31] . It would imply the presence 
of a collapsing cosmic texture at redshift z ~ 6, with 
a typical scale of ^ 5°, corresponding to a symmetry- 
breaking energy scale of ~ 9 x 10^^. In addition, 
some follow-up tests were proposed to confirm or dis- 
card the cosmic texture hypothesis. In particular, one 
should expect to have more spots generated by evolv- 
ing textures in CMB observations (e.g., ~ 7 with scales 
'i^ > 2°); a lack of E-mode polarization is predicted, as 
compared to the values associated with spots derived 
from a Gaussian CMB field; and, finally, a lensing of 
the CMB photons is potentially detectable with future 
small scale CMB experiments as ACT or SPT. 

Summarizing, the study of the Cold Spot has pro- 
vided a wealth of information, and the upcoming high 
quality CMB data, as the one expected from Planck, 
guarantee that still more knowledge will come from the 
analysis of this very interesting feature. 
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